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ABSTRACT 

We investigate statistical distributions of differences in gravitational-lensing deflec- 
tions between two light rays, the so-called lensing excursion angles. A probability 
distribution function of the lensing excursion angles, which plays a key role in esti- 
mates of lensing effects on angular clustering of objects (such as galaxies, QSOs and 
also the cosmic microwave background temperature map), is known to consist of two 
components; a Gaussian core and an exponential tail. We use numerical gravitational- 
lensing experiments in a ACDM cosmology for quantifying these two components. We 
especially focus on the physical processes responsible for generating those two compo- 
nents. We develop a simple empirical model for the exponential tail which allows us 
to explore its origin. We find that the tail is generated by the coherent lensing scatter 
by massive halos with M > 1O 14 /i -1 M0 at z < 1 and that its exponential shape arises 
due to the exponential cut-off of the halo mass function at that mass range. On scales 
larger than 1 arc minute, the tail does not have a practical influence on the lensing 
effects on the angular clustering. Our model predicts that the coherent scatter may 
have non-negligible effects on angular clustering at sub-arcminute scales. 

Key words: gravitational lensing - cosmology: theory - dark matter - large-scale 
structure of universe 



1 INTRODUCTION 

Light rays are deflected when they propagate through an 
inhomogeneous gravitational field, such as the real universe 
we live in. The lensing deflection angle varies from one direc- 
tion to another, and thus the difference in deflection angles 
between two light rays, which we call the "lensing excursion 
angle" , does as well. Consequently, it is not easy to infer the 
transverse distance between two celestial objects at a cos- 
mological distance from their angular separation in the sky. 
Strictly speaking, lacking complete knowledge of the matter 
distribution in the universe, this is impossible to do. 

Since distance is one of most fundamental physical 
quantities, the lack of a precise distance measure to the far 
universe may prevent us from a detailed understanding of 
the universe. A well known case is that angular correlations 
of distant galaxies and of the temperature map of the cosmic 
microwave background are altered by lensing deflections (see 
§9 of Bartelmann & Schneider 2001 for a review and refer- 
ences therein). 

Although we cannot know the lensing excursion angle 
for an individual pair of light rays, knowledge of their sta- 
tistical distribution greatly helps us in estimating the order 



of magnitude of lensing effects. In addition, it allows the in- 
trinsic angular correlation functions to be deconvolved from 
measured correlation functions (Bartelmann & Schneider 
2001). It is thus of fundamental importance to understand 
in detail the statistical distribution of the lensing excursion 
angles. 

The analytic model for computing the variance of lens- 
ing excursion angles in a framework of modern cosmologi- 
cal models was developed by Seljak (1994; 1996), based on 
the linear perturbation theory (the so-called power spectrum 
approach). Hamana & Mellier (2001) performed numerical 
experiments of the gravitational lensing deflections in cold 
dark matter models and examined the statistical properties 
of the lensing excursion angles. They found that the prob- 
ability distribution function (PDF) of the excursion angles 
consists of two components, a Gaussian core and an expo- 
nential tail, and that the variance of the Gaussian core com- 
ponent agrees well with the prediction by the power spec- 
trum approach. They argued that the exponential tail may 
be generated by coherent lensing scattering by massive ha- 
los which is not taken into account in the power spectrum 
approach. 
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The purpose of this paper is twofold: The first is to 
explore the origin of the exponential tail of the lensing ex- 
cursion angle PDF. The second is to develop an empirical 
model for the exponential tail. To pursue these purposes, 
we first examine in detail properties of the exponential tail 
of the excursion-angle PDF using numerical experiments in 
§2. Then in §3, we describe a model for the exponential 
tail which is based on the assumption that the tail origi- 
nates from coherent lensing scattering by individual massive 
dark-matter halos, and compare the model predictions with 
numerical results. We also discuss a general picture of the 
light propagation in the universe paying special attention 
to the role of secular random deflections by either large or 
small-scale structures and a coherent scatter by a massive 
halo. Finally, we give a summary and discussion in §4. 



2 RAY-TRACING SIMULATION 
2.1 VLS iV-body simulation 

We performed weak lensing ray-tracing experiments in a 
Very Large iV-body Simulation (VLS) carried out by the 
Virgo Consortium (Jenkins et al. 2001, and see also Yoshida, 
Sheth & Diaferio 2001 for simulation details). The simula- 
tion was carried out using a parallel P 3 M code (MacFarland 
et al. 1998) with a force softening length of l BO f t ~ 30 /i _1 kpc. 
The simulation employed 512 3 CDM particles in a cubic 
box of 479 h~ Mpc side length, which gives a particle mass 
of m P art = 6.86 x 10 10 1i~ 1 Mq. It uses a flat cosmological 
model with a matter density fio = 0.3, a cosmological con- 
stant J1a = 0.7, and a Hubble constant Hq = 100ft with 
h = 0.7. The initial matter power spectrum was computed 
using CMBFAST (Seljak & Zaldarriaga 1996) assuming a 
baryonic matter density of fib = 0.04. The normalization of 
the power spectrum is taken as <rg = 0.9. 



2.2 Weak lensing ray-tracing simulation 

The multiple-lens plane ray-tracing algorithm we used is de- 
tailed in Hamana & Mellier (2001; see also Bartelmann & 
Schneider 1992 and Jain, Seljak & White 2000; Hamana, 
Martel, Futamase 2000; Vale & White 2003 and Hamana, 
Takada & Yoshida 2004 for the theoretical basics and tech- 
nical issues); thus in the following we describe only aspects 
specific to the VLS simulation data and to ray-tracing ex- 
periments in this study. 

We use thirteen snapshot outputs from two runs of the 
TV-body simulation which differ only in the realization of the 
initial fluctuation field. A stack of these outputs provides the 
density field from z — to z — 6.8. We do not use further 
higher redshift outputs because of two reasons; (1) discrete- 
ness effects of particles (Hamana, Yoshida & Suto 2002), and 
(2) an artificial power excess in the density power spectrum 
due to the "glass" initial condition (White 1996) at around 
the mean separation length of particles, both of them are 
significant at such high redshifts. For higher redshifts up to 
the last scattering surface (z w 1100), we simply consider 
a homogeneous density field. Thus within 6.8 < z < 1100 
rays propagate as in a perfectly homogeneous universe. This 
treatment misses lensing contributions from structures at 



that redshift range. It has turned out that this approxima- 
tion causes only a minor effect. We will discuss its influences 
on our analyses later. 

Each iV-body box is divided into 4 sub-boxes with an 
equal thickness of 119.75 /i _1 Mpc. The A^-body particles in 
each sub-box are projected onto lens planes. In this way, 
the particle distribution between an observer and z = 6.8 is 
projected onto 50 lens planes. Note that, in order to mini- 
mize the difference in redshift between a lens plane and an 
output of iV-body data, only one half of the data (i.e. two 
sub-boxes) of z — output is used. The particle distribution 
on each plane is converted into the surface density field on a 
2048 2 regular grid using the triangular shaped cloud (TSC) 
assignment scheme (Hockney & Eastwood 1988). The grid 
size is 0.23/i _1 Mpc which is chosen to maintain the resolu- 
tion provided by the iV-body simulation and removing at 
the same time the shot noise due to discreteness in the N- 
body simulation (this choice is equivalent to the "large-scale 
smoothing" in Menard et al., 2003, we refer the reader to 
this reference for further examination of the effective resolu- 
tion of the ray-tracing simulation). Its computation follows 
the procedure described in Hamana & Mellier (2001). 

Having produced surface density fields on all lens 
planes, 1024 2 rays are traced backwards from the observer's 
point using the multiple-lens plane algorithm (e.g. Schnei- 
der, Ehlers & Falco 1992). The initial ray directions are set 
on 1024 2 grids with a grid size of 0.25 arcmin, thus the total 
area covered by rays is 4.27 2 square degrees. We produced 
36 realizations of the underlying density field by randomly 
shifting the simulation boxes in the direction perpendicular 
to the line-of-sight using the periodic boundary conditions 
of the iV-body boxes. 

The 36 realizations are not perfectly independent be- 
cause they are generated from the same iV-body outputs 
(but using different combinations of random lines of sight) 
which come from two runs of iV-body simulation. Therefore 
the generated lensing data (the lensing deflection field, lens- 
ing convergence and shear map) are subject to sample vari- 
ance. In order to test its magnitude, we compare the conver- 
gence two-point correlation function with its theoretical pre- 
diction (Jain & Seljak 1997) in Figure^ The measurements 
from the ray-tracing experiment are plotted by symbols with 
error bars which represent the mean and root-mean-square 
among the 36 realizations, while the solid lines show the 
prediction. Note that the measurement for z s = 1100 should 
be compared with the dotted line which shows the theoret- 
ical prediction for z s = 1100 but the contribution from the 
density fluctuations between z = 6.8 and 1100 is ignored. 
The measurements are in good agreement with the theoreti- 
cal prediction in shape, but are slightly higher in amplitude. 
This excess may be mostly attribute to the sample variance 
and implies that there exists an excess power in the lensing 
potential field. Since lensing deflections result from the same 
potential field, it is expected that there exists, to a similar 
extent, an excess in the deflection angle statistics. On scales 
smaller than 1 arcmin, the slope of the measured correlation 
function becomes flatter than predicted; this is due to the 
limited resolution of the TV-body simulation. The effective 
angular resolution of the convergence field is about 1 arcmin 
for lower redshift (z s < 3) and is slightly better for higher 
redshifts (see Menard et al., 2003 for further discussion on 
the resolution issue). 
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Figure 1. The convergence two-point correlation functions. Mea- 
surements from ray-tracing experiments are shown by symbols 
with error bars which represent the mean and root-mean-squarc 
among the 36 realizations. Crosses, filled triangles and open circle 
are for z s = 1, 3 and 1100, respectively. The error bars of the cases 
for z a = 3 and 1100 are not displayed for clarity but are as similar 
magnitude as the z B = 1 case. The solid lines show the theoretical 
prediction (e.g., Jain & Seljak 1997) in which the fitting function 
of nonlinear power spectrum by Peacock & Dodds (1996) is used 
to include the effect of the nonlinear growth of the density field. 
The dotted line shows the theoretical prediction for z s = 1100 
but the contribution from density fluctuation at 6.8 < z < 1100 
is not integrated. 



2.3 PDF of the lensing excursion angles 

Using weak-lensing experiments, we study the statistics of 
differences in deflection angles between two light rays, which 
we refer to as the "lensing excursion angle". The deflection 
angle of a light ray, which is computed by the lens equa- 
tion, is simply the difference between its positions 8i and 
8s on the image and source planes, respectively. Denoting 
the deflection angle of two rays by a 1 and a 2 , respectively, 
we write the lensing excursion angle between these rays as 
5a — a 1 — a 2 . Similarly we denote their intrinsic separation 

b y e 1 2 = \e 1 s -o%\. 

Let us first look into the probability distribution func- 
tion (PDF) of the lensing excursion angles which is one of 
most fundamental statistics. Figure|2|shows the PDFs of the 
lensing excursion angles normalized by its intrinsic separa- 
tion (i.e. 5a/ 812). Since the vector field 5a has no preferred 
direction, we use both components, Sai and 5ct2, to com- 
pute the PDFs. The dotted lines in each plot of Figure [5] 
show the Gaussian PDF with its standard dispersion (a 2 ) 
computed from the PDF itself (i.e. a 2 = J dx i 2 PDF(a;)). 
As was first pointed out by Hamana & Mellier (2001), the 
PDFs consist of two components, a Gaussian core and the 
exponential tail, which are generated by different physical 
processes as we will discuss below. 

The origin of the Gaussian core is explained as follows: 
Light rays from a cosmological distance undergo many (ei- 
ther strong or weak) gravitational lensing deflections. Since 
the spatial distribution of lenses at a large separation are un- 
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Figure 2. The solid curves show the probability distribution 
function of the lensing excursion angles normalized by the intrin- 
sic ray separation. The source redshift is z s = 1100. The ranges of 
the intrinsic separation of the light-ray pairs 812 are given in each 
panel. The dotted curves show Gaussian distributions with their a 
computed from the measured PDFs (i.e., a 2 = J dx i 2 PDF(i)). 



correlated, rays basically undergo many uncorrelated deflec- 
tions. Provided the separation between two rays is so large 
that effects of coherent scattering can be ignored, two light 
rays undergo independent deflections. According to the cen- 
tral limit theorem, the statistical distribution of the lensing 
excursion angles of such light-ray pairs is given by a Gaus- 
sian. A necessary condition for the central limit theorem to 
hold is that the parent distribution of the individual events 
which are being superposed has finite variance. The deflec- 
tion angle calculated in the weak-lensing regime using the 
power-spectrum approach has finite variance (Seljak 1994; 
1996), but it is based on linearized gravity and ignores lens- 
ing by individual halos. On the other hand, numerical grav- 
itational lensing experiments show that the PDF consists of 
the Gaussian core and the exponential tail. In order to un- 
derstand whether the Gaussian core can indeed be caused 
by the superposition of many deflections, we need to inves- 
tigate the variance of the excursion angle. In particular, the 
numerical experiments miss the influence of numerous dis- 
tant lenses, because of their necessarily finite volume. We 
will now show, using a simple approach, that accumulated 
contributions from very distant lenses do not significantly 
affect the excursion angle variance, thus it remains finite. 

Consider first a single light ray passing the lens plane 
in the origin. There is a finite number of lenses close to the 
ray, thus we can restrict ourselves to distant lenses since 
we are investigating whether the deflection-angle variance 
is finite or not. Axially symmetric lenses more distant than 
their (e.g. virial) radii act as point lenses, thus we can ap- 
proximate their deflection angles by oti = 8i/8f. Assum- 
ing the lenses have a number density n and are randomly 
distributed, the variance of the total deflection angle con- 
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Figure 3. The standard dispersion of the lensing excursion an- 
gles plotted as the function of the separation angle. Symbols and 
error bars show the mean and RMS among 36 realizations of the 
ray-tracing numerical experiments. Solid lines represent the the- 
oretical predictions from the power spectrum approach (Seljak 
1996), in which the fitting function of the nonlinear power spec- 
trum by Peacock &; Dodds (1996) was used. The source redshifts 
are, from lower to higher, z B =l, 3, 6.8 and 1100, respectively. 
The dotted line shows the theoretical prediction for z B = 1100 
but the contribution from density fluctuation at 6.8 < z < 1100 
is ignored. 



tributed by lenses in a ring around the origin with radius 9 
and width d9 is 



2im0 d9 



2nn d In f 



(1) 



Integrating over 9 shows that the variance diverges logarith- 
mically. 

The situation changes for the excursion angle. Con- 
sider two light rays piercing the lens plane at positions 
01,2 = (=Fd/2,0). Specializing again to distant lenses, we 
can approximate the individual deflection angles by those of 
point lenses. The excursion angle of the lenses in a ring of ra- 
dius 6 and width d8 around the origin can then be expanded 
to lowest order in d/9, 



N d 

Sa(9) = 22 ~{p (-cos2<^,sin2<^) , 

i=l 



(2) 



where fa is the polar angle of the i-th lens. Again assuming 
randomly distributed lenses, the variance of the excursion 
angle contributed by the lenses in the ring is thus 



(Sa 2 (9)) = 




2n 2 d 2 n9 d9 



(3) 



i.e. the excursion-angle variance converges like 9 when 
integrated over 9 to infinity. Thus, we can apply the central 
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Figure 4. The PDF of lensing excursion angles normalized by 
its intrinsic separation (plotted positive side only). The source 
redshifts are, from narrower to broader PDF, z s =X, 3 and 1100, 
respectively. As these plots show, for a given 9 range, the ampli- 
tude of the PDF tail becomes higher as z s increases, while their 
slope is almost unchanged for plotted redshifts z s > 1. 



limit theorem to the excursion angle, while we could not for 
the deflection angle itself. 

Now we test the theoretical model prediction for the 
variance of the excursion angles developed by Seljak (1994; 
1996) against our numerical results. Figure [3] compares the 
standard dispersion measured form the numerical experi- 
ments with the theoretical prediction. The dispersion be- 
comes larger as the light rays travel a longer distance, be- 
cause the rays can undergo more deflections. It is found in 
the plot that the measurements are slightly larger than the 
prediction. However, a similar excess is seen in the conver- 
gence correlation function (Figure , thus this is mostly 
due to the sample variance. We may therefore conclude that 
the power spectrum approach provides a good prediction 
even for z s = 1100, and the non-Gaussian tail has no strong 
contribution to the variance. It is important to notice that 
coherent scattering by lensing due to massive halos that gen- 
erate the exponential tail contribute only very little to the 
excess in the measured dispersion over the prediction. 

Let us now turn to the tail of the lensing excursion 
angle PDF. Figure |3] compares the PDF obtained from the 
ray-tracing numerical experiments for three source redshifts, 
z s = 1, 3 and 1100, and for various ranges of ray separa- 
tions. This Figure represents major characteristics of the 
non-Gaussian tail: (a) it has an approximately exponential 
slope; (b) it changes little with the source redshift, but its 
amplitude increases with the source redshift, at least within 
the redshift range we consider [z 3 > 1). We fit the tail of 
the PDFs to the exponential distribution: 



E(x) = pexp(-qx). 



(4) 



To do this, we take two points x\ and X2 such that 



PDF(xl) = 1.0 x 10" 2 and PDF(x2) 



1.0 x 10" 



The 
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Figure 5. Parameters in the exponential distribution (eq.[IJ com- 
puted by fitting to the non-Gaussian tail of PDFs from ray-tracing 
experiments. The filled circles and crosses show the results from 
the z s = 1100 and z s = 3 cases, respectively. The error bars rep- 
resent the ranges of ray separations taken to compute the PDFs. 



comparison between the results of z s = 3 and z s = 1100 
plotted in Figure confirms the above point (b) in a quan- 
titative manner. Note that fitting the exponential function 
to the PDF tails becomes poor for large ray separations be- 
cause the non-Gaussian tail does not appear prominently 
due to the limited statistics. This accounts for the steep rise 
of both p and q at large 612 which rather reflects the slope 
of the Gaussian core. 

In the next section, we develop a model of the expo- 
nential tail and explore its origin. To do this, a visual im- 
pression from Figure |n] could be informative. In this figure, 
the celestial distributions of ray pairs having a large lensing 
excursion angle and of massive halos (Mhaio > 10 13 /i _1 M©) 
in one realization of the numerical experiment are displayed 
(see Hamana et al. 2004 for a detailed description of the 
construction of a halo catalog on a light-cone). Only the ha- 
los within the redshift interval between and 1 are plotted 
using black symbols. Red dots represent the ray pairs which 
obey the following criteria, an unlensed ray separation of 
#12 > 1 arcmin and an excursion angle of \Sct\/9i2 > 0.7. 
Apparently, most of large excursion angle ray pairs pass very 
close to a massive halo. We argue in the following section 
that the exponential tail results from coherent strong deflec- 
tions of two nearby rays by a massive halo, and explain the 
origin of the above characteristics using simple models. 



3 ORIGIN OF THE EXPONENTIAL TAIL OF 
THE EXCURSION ANGLE PDF 

In this section we explore the origin of the exponential tail 
of the excursion angle PDF found in the ray-tracing nu- 
merical experiments. For this purpose, we focus on the tail 




Figure 6. The celestial distributions of ray pairs having a large 
lensing excursion angle and of massive halos in one realization of 
the numerical experiment. See Hamana et al. (2004) for a detailed 
description of construction of the halo catalog on a light-cone. 
Black symbols represent halos; the large filled circles, small filled 
circles and dots are for halos with M halo > 4 X 10 14 , 4 X 10 14 > 
A4aio > 4 x 10 13 and 4 x 10 13 > M halo > 1 x 10 13 [ft- 1 Af ], 
respectively. Only the halos within the redshift interval between 
and 1 are displayed. Red dots represent the ray pairs which 
obey the following criteria, an unlensed ray separation of 812 > 1 
arcmin and an excursion angle of \6oi\/6v2 > 0.7. The middle 
points of ray pairs are displayed. 



part and do not consider the Gaussian core whose origin 
has been investigated in the literature (Seljak 1994; 1996; 
see also chapter 9 of Bartelmann & Schneider) and also in 
the last section. We develop a theoretical model from two 
assumptions: Large excursion angles are mainly caused by 
the strong lensing of a massive halo, and the probability for 
a ray to undergo multiple strong lensing events is negligi- 
ble. The former is reasonable because a process that is not 
taken into account in the power spectrum approach could 
generate non-Gaussian features. Also the visual impression 
from Figure [(J could be a support of that idea. The latter is 
validated by the observational fact of the small cross section 
for strong lensing events by a single lens (either a galaxy 
or cluster of galaxies) such as multiply-imaged QSOs and 
strongly-lensed arc-like images of distant galaxies. Thus, it 
is certain that multiple scattering by more than one massive 
halo is very rare. 

We consider the same ACDM cosmology as one adopted 
for the numerical experiments in §2. We denote the PDF for 
finding a ray pair with 812 having the excursion angle Set 
by PDF(<5a|0i 2 ). 

3.1 Lensing deflection by a universal density 
profile halo 

Navarro Frenk & White (1996; 1997, NFW hereafter) found 
from A r -body simulations that the density profile of dark 
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matter halos can be fitted by a universal form regardless 
of their mass and redshift. We adopt a truncated universal 
profile; 



p(x) 



x s (l + z) 3 



v 



(•») 



for r < rvi r and otherwise, where r s and r v j r are the scale 
radius and virial radius, respectively. It is convenient to in- 
troduce the concentration parameter c v ir = r v \i/r s . Navarro 
et al. (1996) proposed the universal inner slope of s = 1, 
while a steeper slope was claimed by later studies using 
higher resolution iV-body simulations; Moore et al. (1998; 
1999), Ghigna et al. (2000) and Fukushige & Makino (2001; 
2003) found larger values such as s — 1.5, while Jing (2000) 
and Jing & Suto (2000) pointed out that it varies from 1.1 to 
1.5 and argued a possible weak dependence on the halo mass. 
In this paper, we consider two cases s — 1 and 1.5. Navarro 
et al. (1997) and Bullock et al. (2001) have extensively ex- 
amined a relation between the concentration parameter and 
the halo mass and its redshift evolution adopting a fixed 
value of s = 1. We adopt a generalized mass-concentration 
relation proposed by Oguri, Taruya & Suto (2001; see also 
Keeton & Madau 2001); 



r (M,z) = {2-s)- 



M 



1 + 2 V W^h^Mp 



(6) 



Bullock et al. (2001) suggested c, ~ 8 for the ACDM model, 
which we adopt as a fiducial choice. We note that there is a 
relatively large scatter in this relation (Bullock et al. 2001; 
Jing 2000). The virial mass (defined by the mass within the 
virial radius r v i r ) of the universal halo is given by 



•? m(c v ii,s) 

r" 3 



with 

m(Cv ir , s) = 



dx 



(l + xf- 



(8) 



Since the spherical collapse model indicates that M v \i = 
47rrvi r <5 v ir(2)po/3, where 8 V1T is the over-density of collapse 
(sec Nakamura & Suto 1997 and Henry 2000 for useful fitting 
functions), one can express p a in terms of 8 v ir(z), c v i r and s: 



Ps = 



r. _ 3 

OvirPO C vir 

3 m(c v ir,s) 



(9) 



Let us summarize basic equations for gravitational lens- 
ing properties of the truncated universal profile halo (Takada 
& Jain 2003, see Bartelmann 1996; Wright & Brainerd 2000; 
Oguri et al. 2001 for lensing properties of the non-truncated 
universal profile lens model). The surface mass density of 
the truncated universal profile halo is given by, 



" \/ C vir — y2 

£(!/)=/ Z dz p(y,z) = 2p s r s f(y), y=-, (10) 



with 



1 



f(y) = / dz 

Jo (y 2 + z 2 y/ 2 (l + vV + z 2 )^) 



(11) 



for y ^ c v i r and f(y) — otherwise. The projected mass 
within a radius b is 

M(< b) = 2n f dy' y'E(y') 



4nr a p a / dy yf(y), 



(12) 



where h — b/r s . We perform the above integration numeri- 
cally. The thin lens equation is written by 



6s = 6i - a{Di0i), 
with 

... 4GM(< b) D U 
a{b) = 



(13) 



(14) 



c 2 b D s ' 

In the last expression, the origin of the coordinates is taken 
at the lens center, b = Didi is the impact parameter, and 
Di, Dis and D s are the angular diameter distances from 
observer to lens, from lens to source, and from observer to 
source, respectively. The deflection angle of the truncated 
universal profile halo is given by, 



a(0) = a*g(e), 
with 



and 



Jo dy yf(y) 



D,6 



m(Cvir,s) X 



(15) 



(16) 



20, 



Dis ( Ho 



(Hoy 2 , 



^vir 



lft-iMpc I D 



(17) 



Note that a* oc M 2 ^ 3 . It is important to notice that a de- 
pendence of the halo profile parameters on the deflection 
angle enters only through the function g{9). Note that for 
6 ^ # v ir it reduces to g(9) — 8 v i r /8 (where # vlr is the angu- 
lar virial radius defined by (9 v ir = r v i r /Di). In Figure |7| the 
function g(0) is plotted for various value of c*. As one may 
see in the Figure, the deflection angle profile g{6) peaks at 
9 ~ ^ v i r /c, = 8 S (0 B = r s /D{) and the peak value does not 
strongly depend on the inner slope s. It is also found that the 
peak value relates to the concentration parameter roughly 
by Smax ~ 0.1c» + 1. Therefore, in a reasonable range of 
c* the maximum deflection angle by a single universal halo 
lens is a max = (1 — 3)a„. One may also find that at the 
inner part, the deflection angle is larger for a larger c* or for 
a steeper inner slope (thus for more centrally concentrated 
halos). We find that g(0) has an asymptotic inner slope of 



(tx 



*) for s = 1 (s = 1.5). 



3.2 Lensing excursion angles 

Since the deflection angle of a universal density profile halo 
is finite, the excursion angle is finite as well. Clearly, the 
largest excursion angle is 2a max which happens when one 
ray passes at the distance ~ 9 S from the lens and the other 
ray passes at the same distance in the opposite side of the 
lens. Thus this happens only if 812 = 28 s , and is very rare. 
Let us consider the maximum excursion angle produced for 
other ray separations. If 6*12 > 8 S , the largest excursion angle 
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Figure 7. Deflection angle profiles of the truncated universal 
density profile lens, g(9) defined in eq. 1161 . as a function of the 
impact parameter (normalized by the angular virial radius) for 
various concentration parameters denoted in the plot. The left 
panel is for s = 1 and right panel for s = 1.5. Note that g(6) has 
an asymptotic inner slope of cx 0°- 82 (oc - 48 ) for s = 1 (s = 1.5). 



is in the range a max < Sa < 2 a max . This happens when one 
ray passes at ~ 9 S and the another ray passes at the opposite 
side of the lens in the direction connecting the lens center 
and the first ray. While if O12 < 6 S , the largest excursion 
angle is smaller than a max . An important consequence of 
this is that for ray pairs with the separation angle larger 
than 8 3 , the maximum excursion angle does not strongly 
depend on the ray separation but lies in a small range of 
(1 — 2)a max , and it scales with the halo mass as oc M 2 ^ 3 . 



3.3 PDF of Sa 

Let us first consider the probability distribution induced by 
one halo, which we denote by PDFi(<5a:|#i2). Let p(6) be 
the probability of a ray passing at a small area — > + 60 
from a lens center, which is given by the cross section area 
(denoted by A) normalized by the unit solid angle (dQ,): 



P(0) 



A{6 -^0 + 56) 



(18) 



Then PDFi((5a|6'i2) is given by the joint probability, 
d 2 t 



PDFi(«5a|0i 2 ) 



2-7T 



-p(0i) p(0i +01 2 ), (19) 



where 0i2 = {#12 cos(0i2), #12 sin(<^i2)}. Note that since we 
are considering lensing by a single halo having a certain 
density profile and mass, given a configuration of a light ray 
pair, its excursion angle is uniquely determined. The total 
PDF((5a|#i2) is obtained by summing PDFi (5a|#i2) over 
halos within a light-cone volume, 

PDF{Sa\9 12 ) = J dV J dM n ha i (M, 2 )PDFi(5a|<9i 2 ) 

r T ( z s) 

dr r 2 I dM 



10- 
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xn halo (M,r[«])PDFi(<5a|<9i 2 ), (20) 



Figure 8. PDF of the lensing excursion angles from one halo, 
PDFi(<5a), for three halos masses M halo = 10 13 , 10 14 and 
1O 15 /i _1 M0. The normalization is arbitrary. The lens and source 
redshifts are z; = 0.3 and z s = 1100, respectively. The solid lines 
are for s = 1 and dashed lines are for s = 1.5. The concentration 
parameter taken is c* = 8 for all cases. 

where r is the comoving radial distance, dV = r 2 dr (this 
expression is valid only for a fiat cosmological model) is 
the unit volume element and nh a i (M, z) is the halo mass 
function. We adopted the mass function by Sheth & Tor- 
men (1999). Note that this approach breaks down for small 
excursion angles where secular scattering by distant and/or 
small halos are important, which generate the Gaussian core 
of the excursion-angle distribution. 



3.4 Results 

Let us start with the excursion angle PDFs from one halo 
plotted in Figure |H| which help to understand the origin of 
the exponential tail. The most important point which should 
be noticed is the sharp cut-off in a large excursion angle. 
This is a natural consequence of the fact that the deflection 
angle of the universal density profile halo is finite (see 
The maximum excursion angle scales with the halo mass 
roughly as oc M 2 ' 3 (with a small correction by the mass 
dependence of the concentration parameter) as far as the 
separation angle is larger than the angular scale radius (9 S ) 
of a lensing halo, as explained in the last subsection. The 
other important point is that the mass-independent double 
power-law slope, its power-law slope is oc 5oT for smaller 
excursion angles and oc <5a -2 ' 64 for larger angles. The former 
is generated by ray pairs in which both rays pass outside of 
the virial radius, while the latter is generated by pairs of 
rays of which one passes outside of the halo, and the other 
inside. 

Under the assumptions stated in the last subsection, 
the excursion angle PDF is obtained by summing up con- 
tributions from single halos over a wide range of the halo 
mass and integrating over the redshift of halos as defined 
by eq. 1201 . We plot the PDF computed from such a model 
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Figure 9. Model prediction of the lensing excursion angle PDF 
for z s = 1100, c* = 8 and s = 1. The black line shows the full 
PDF, while the colored lines show contribution from a limited 
range of the halo mass denoted in the panel. 
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Figure 10. Model prediction of the lensing excursion angle PDF 
for z s = 1100, c* = 8 and s = 1. The black line shows the full 
PDF, while the colored lines show contribution from a limited 
range of the lens redshift denoted in the panel. 



in Figure |2l The black line shows the total PDF, while col- 
ored lines represent contributions from narrow limited mass 
ranges. This Figure clearly illustrates the origin of the ex- 
ponential tail. There are two key points; a large excursion 
angle can only be generated by massive halos with mass 
typically larger than 10 14 h~ 1 Mo, and at such a halo mass 
range, the mass function decreases exponentially. Accord- 
ingly, the number of more massive halos that can contribute 
to a larger excursion angle decreases exponentially, and as 
a result, the exponential slope of the PDF arises. 

In order to examine what redshift range of halos makes a 
major contribution to the exponential tail, we plot in Figure 
I1UI the excursion angle PDFs computed for limited ranges 
of the lens redshifts (upper panel) and their percentage of 
the contribution (lower panel). It is seen in the lower panel 
that contributions from halos at redshifts below 1 account 
for almost full amplitude of the PDF tail. This is explained 
by a rapid evolution of the halo mass function at the high 
mass end. In fact, the number density of massive halos with 
A^haio ^ 2x 10 14 ft _1 M0 decreases by more than one order 
of magnitude from z = to 1. 

The above results lead to the following explanation for 
the origin of the lensing excursion-angle PDF which has two 
components, a Gaussian core and an exponential tail: The 
light rays emitted at high redshifts undergo many gravi- 
tational deflections by (large- or small-scale) structures on 
their way to us, which after many uncorrelated deflections 
produces the Gaussian core. Some small part of the rays 
is strongly lensed by a massive halo with mass larger than 
10 14 h~ x M© at a low redshift of 2 < 1, and the coherent 
deflections caused by the strong lensing produce the expo- 
nential tail. Therefore even if a ray pair encounters a strong 
coherent deflection by a single massive halo, its excursion 
angle is not solely determined by the strong halo lensing 



but the random scattering also contributes to it to a smaller 
extent. This effect on the PDF is taken into account by the 
convolution: 



E{x)= / dyG{y)E'{x-y), 



(21) 



where G(y) denotes the Gaussian distribution, and E'(x) 
denotes the tail part produced by the halo lensing with- 
out considering the contribution from the random scatter- 
ing (which could be computed by the model described in 
this section) . Since, as shown in Figure |3J the tail is well 
approximated by the exponential shape with a constant in- 
dex over a wide range of 5a, it is reasonable to approximate 
E'(x) ~ p'exp(— qx) and, 



.»„ r / dy ex P[-7^5 I exp(-q(x - y)) 



E(x) 

= pexp(-ga;), 
with the boosted amplitude, 

p — p exp ' 



(22) 



(23) 



Thus under the above approximation, the slope of the ex- 
ponential tail is unchanged but the amplitude is increased. 
This combined with the fact that almost all contributions 
to the PDF tail come from coherent scattering by massive 
halos at z < 1 (Figure accounts for the trend observed 
in the excursion angle PDFs obtained from the ray-tracing 
simulation that its slope does not depend strongly on the 
source redshift but its amplitude increases with the source 
redshift. Actually, a very similar trend is observed in the 
model PDFs plotted in Figure lTTI which shows the corrected 
PDF tails for z 3 = 1, 3 and 1100 and for four ray separations. 
Here, in order to compute the boost factor of eq. 12MI) . we 
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Figure 11. The model predictions of the lensing excursion angle 
PDF after the correction to the accumulative lensing effect being 
made (using the approximate way of eq. The halo param- 

eters are c* = 8 and s = 1. The ray separations are denoted in 
each plot. The source redshifts are, from inner to outer, z 3 =l, 3 
and 1100, respectively. 



compute p' and q by fitting the model PDFs at two points, 
PDF(kI) = 1.0 x 1(T 2 and PDF(a;2) = 1.0 x 10~ 3 to the 
exponential function. Note that after this correction the am- 
plitude of the tails can be increased more than one order of 
magnitude because of the steep slope of the exponential tail 
(thus for a large g). 

Four panels of Figure 1121 show the corrected model 
PDFs for various values of c* and s, which help to under- 
stand the dependences of the shape of the PDF tail on the 
halo parameters. Clearly, the broader tail appears for models 
with a larger c, or a larger s, because such models generate 
a larger maximum deflection angle. It is important to notice 
that a small change of c* or s causes a very large, nonlinear 
change in the shape and amplitude of the PDF. Therefore, 
choosing mean values of c* and s does not provide a mean 
PDF but gives a lower amplitude. It should be noted that 
the halo model parameters indeed have a large scatter (Jing 
2000; Jing & Suto 2000; see also Figure 9 of Hamana et 
al. (2004) which clearly shows that the compactness of halo 
mass distributions has a large scatter). 

Finally, we compare in Figure 1131 the parameters in 
the exponential function eq. 1122ft measured from the model 
PDFs (symbols with lines), with the results from the ray- 
tracing experiments. We evaluate p and q by fitting the cor- 
rected model PDF (i.e., after the correction by eq. Il22t be- 
ing made) . As shown in Figure 1131 the parameter q corre- 
lates with the slope parameter q as expected. The measured 
exponential slope parameter q plotted in the lower panel 
are larger than the results from the numerical experiments, 
though the slope with the separation angle is very similar. 
The discrepancy is smaller for models with a larger c* or a 
larger s. Therefore the model prediction may be improved 








Figure 12. The model predictions of the lensing excursion angle 
PDFs. We have applied the corrections given by equations 1221 
and 1231 . Ray separations are 9 = 2, 8 and 32 arcmin (solid, 
dotted and dashed, respectively). The source redshift is taken by 
z s = 1100. Halo model parameters are denoted in each panel. 
The thin lines in the top right and two bottom panels show, for 
comparison, the predictions of the fiducial model plotted in the 
top left panel. 



if one takes into account the scatter in the halo parameters 
c» or a larger s. 

Quantitatively, none of the four models plotted in Fig- 
ures ^| and ED are i n very good agreement with the simu- 
lation results. The discrepancy is partly due to the scatter 
in the halo model parameters as has been discussed above. 
Also, a deviation in the halo mass distribution from spheri- 
cal symmetry could partly account for it. Actually, the mass 
distribution of most of the halos significantly deviates from 
spherical symmetry (Hamana et al. 2004). In addition, the 
spatial correlation of halos may have an influence on the 
excursion angle PDF, because massive halos are strongly 
clustered. A close look at the sky distributions of ray pairs 
having a large excursions angle and of massive halos shown 
in Figure |H| reveals that those deviations from our simple 
model should indeed play a role; namely, it is seen in the 
Figure that a small part of most massive halos does not 
produce a large lensing excursion event, and that a small 
part of large excursion angle ray pairs does not intersect a 
very massive halo. 

We may conclude, from what has been seen above, that 
our simple model succeeds in getting the essential mecha- 
nism of generating the exponential tail and in explaining 
the origin of the major characteristics of the tail. The model 
predictions are in reasonable agreement with the simulation 
results. Further modifications of the model taking into ac- 
count details of halo properties, such as scatter in the halo 
model parameters, deviations of the halo mass distribution 
from spherical symmetry and clustering of halos, are needed 
to improve the accuracy of the model prediction.. 
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Figure 13. Parameters in the exponential distribution (eq. 1221 
computed from the fit to the model predictions of excursion angle 
PDFs for z s = 1100. Model parameters are denoted in the plot. 
Dotted lines show the results from ray-tracing simulation. 



4 SUMMARY AND DISCUSSIONS 

We have investigated the statistical distribution of lensing 
excursion angles, paying a special attention to the physical 
processes that are responsible for generating two compo- 
nents of the PDFs: the Gaussian core and the exponential 
tail. We have used the numerical gravitational lensing ex- 
periments in a CDM cosmology to quantify these two com- 
ponents. 

The origin of the Gaussian core is explained by the ran- 
dom lensing deflections by either linear or nonlinear struc- 
tures. The variances of the Gaussian core measured from 
the results of the numerical experiments are found to be in 
a good agreement with the prediction by the power spectrum 
approach (Seljak 1994; 1996). 

The presence of the exponential tail was first found 
by Hamana & Mellier (2001) but its origin remains unre- 
vealed. The tail is characterized by two parameters: the 
slope and amplitude. We have found from the numerical 
experiments that the slope changes little with the source 
redshift while the amplitude becomes greater as the source 
redshift increases, at least within the redshift range we con- 
sider (1 < z 3 < 1100). Since the random lensing deflections 
result in the Gaussian core, the exponential tail is most likely 
to result from coherent deflections. In addition, in order to 
generate a large excursion angle, massive virialized objects 
should be responsible for the exponential tail. Therefore, we 
supposed that the exponential tail originates from coherent 
lensing scatters by single massive halos. 

We have developed a simple empirical model for the 
exponential tail of the lensing excursion angles PDF. We 
used the analytic models of the dark matter halos, namely 
the modified Press-Schechter (Press & Schechter 1974) mass 



function (Sheth & Tormen 1999) and the Universal density 
profile first proposed by Navarro et al. (1996). Although we 
only consider a coherent lensing scatter by a single massive 
halo and did not take into account the scatters in the halo 
parameters (the concentration parameter c, and the inner 
slope s), our model reasonably reproduces the exponential 
tails computed from the numerical experiments. It is found 
that the massive halos with M > 10 14 /i _1 Mq are responsi- 
ble to the tail and that the exponential shape arises as a con- 
sequence of the exponential cutoff of the halo mass function 
at such mass range. Almost all contributions to the tail come 
from the halos at redshifts below 1. Therefore, the slope of 
the tail is formed by the halos at z < 1. On the other hand, 
the amplitude of the tail is determined by the convolution 
of two contributions, the coherent scatter and the random 
deflections. Since the contribution from the random deflec- 
tions becomes greater as the source redshift increases, the 
amplitude of the tail becomes greater for a higher source 
redshift. These explain the redshift-independent slope and 
the redshift-dependent amplitude found from the numerical 
experiments. 

Does the exponential tail have an influence on the an- 
gular power spectrum of the temperature map of the cos- 
mic microwave background (Ci) ? As far as angular scales 
larger than 1 arcmin are concerned, the answer is no. When 
one computes the lensed Ct one can safely use the approx- 
imate convolution equation given by Seljak (1996), because 
the key assumption in the approximation made for deriving 
the convolution equation is not the Gaussianity of the ex- 
cursion angle PDF, but that its variance is small (cf. §9 of 
Bartelmann & Schneider 2001). Since the smallness of the 
variance is also the case for lower redshifts, the same convo- 
lution technique can be applied to other angular correlation 
functions such as that of galaxies and QSOs. 

Before closing this paper, we present predictions for tails 
of the excursion angle PDF for sub-arcmin ray separations. 
We plot in Figure PHI our empirical model predictions. It is 
found that for ray separations #12 > 10 arcsec, the ampli- 
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tude of the tail keeps increasing with the slope becoming 
flatter in a similar rate of larger separations. However below 
that separation, the growth rate becomes smaller gradually. 
The standard dispersions of these distributions are smaller 
than unity but can be of order O(0. 1) which is comparable to 
that of the Gaussian core, though our simple model adopting 
average halo parameters tends to predict greater amplitude 
than the results from the numerical experiments (see §3.4). 
Therefore, it is possible that on arcsecond scales coherent 
lensing deflections have non-negligible influence on the an- 
gular clustering of objects in the distant universe. Note that 
even if taking the exponential tail into account, the stan- 
dard dispersion of the excursion angles is less than unity, 
thus the approximate convolution equation can be still valid 
but the contribution from the tail to the dispersion should 
be included. We notice that it is however not clear whether 
the assumptions in our model are still valid on such small 
ray separations. The statistical distribution of lensing ex- 
cursion angles for arcsecond separation ray pairs should be 
investigated in a future work with a gravitational numerical 
experiment having a higher resolution. 
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